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ABSTRACT 

Recent observations show evidence that high-z (z ^ 2 — 3) early-type galaxies (ETGs) 
are more compact than those with comparable mass at z 0. Such a size evolution is 
most likely explained by the 'Dry Merger Sceanario'. However, previous studies based 
on this scenario are not able to consistantly explain both the properties of the high- 
z compact massive ETGs and the local ETGs. We investigate the effect of multiple 
(~| , sequential dry minor mergers on the size evolution of the compact massive ETGs. 

^ f From an analysis of the Millennium Simulation Database, we show that such minor 

(stellar mass ratio M2/M1 < 1/4) mergers are extremely common during hierarchical 
structure formation. We perform N-body simulations of sequential minor mergers with 
parabolic and head-on orbits, including a dark matter component and a stellar compo- 
nent. Typical mass ratios of the minor mergers arc 1/20 < M2/M1 ^ 1/10. We show 
that sequential minor mergers of compact satellite galaxies are the most efficient at 
promoting size growth and decreasing the velocity dispersion of the compact massive 
ETGs in our simulations. The change of stellar size and density of the merger rem- 
nants is consistent with recent observations. Furthermore, we construct the merger 
histories of candidates for the high-z compact massive ETGs using the Millennium 
Simulation Database, and estimate the size growth of the galaxies by the dry minor 
merger scenario. We can reproduce the mean size growth factor between z = 2 and 
z = 0, assuming the most efficient size growth obtained during sequential minor merg- 
ers in our simulations. However, we note that our numerical result is only valid for 
merger histories with typical mass ratios between 1/20 and 1/10 with parabolic and 
head-on orbits, and that our most efficient size growth efficiency is likely to an upper 
limit. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxy: 
formation - galaxies: kinematics and dynamics - galaxies: structure - methods: A^- 
body simulations 



1 INTRODUCTION 

There is significant observational evidence that massive 
(> 10^^ Mq) early-type galaxies (hereafter ETGs) at high-z 
(2: ~ 2 — 3) are more compact than galaxies with comparable 
mass in the local Universe. However, a physical explanation 
of the size evolution from such compact ETGs to the local 
ETGs is still lacking within the ACDM galaxy formation 
framework. 

Many observations of the high-z ETGs show that the 
average size of the high-z massive ETGs is smaller by a fac- 
tor of three to five than that of the local ETGs with com- 
parable mass fe.g. iTruiillo et all I2OO6I. iTruiillo et al.ll2007l. 



iBuitrago et al.il2008l . lCimatti et al.ll2008l . lvan Dokkum et all 
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l2008h . These high-z ETGs also lie below the stellar mass-size 
relation in the local Universe; a relation well investigated 
through the Sloa n Digital Sky Survey (SDSS) database (e.g. 
IShen et al.ll2003l 'l. Additionally, observed high-z ETGs indi- 
cate that stellar densities of the high-z ETGs are 1 — 2 or- 
der of magnitude higher than the local ETGs. Several stud- 
ies claim that this size gr owth is even more sig nificant for 
more massive ETGs (e.g. iNewman et aLllioio i. Recently, 
several observations have measured the velocity disper- 
sions of the high-z ETGs (e.g. Ivan Dokkum. Kriek. fc Franxl 
liooi). The observed velocity dispersions are also higher 
than the average velocity dispersion of the local ETGs with 
comparable mass (ICenarro fc Truiilldl2009l.lCappellari et al.l 



l200g . Ivan Dokkum et al.l |2009| . Ivan de Sande et al.1 bOllh 

The compact massive ETGs at high-z are already quiescent 
in star formation activity and have an old stellar population 
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fe.g. lCimatti et al.|[2008l , lvan Dokkum et alfeOOSl) sim ilar to 
ETGs in the local Universe (e.g. Ixhomas et al.lbOOsI '). 

Such comp act massive ETGs are extremely rare in the 
local Universe. iTavlor et al.1 (|2010l ) have searched for com- 
pact ETGs in the local Universe with the SDSS database 
but have only found less massive compact ETGs which 
ar e not readily c ompa rable to the ETGs found at high- 
z. ICimatti et alj (|2008l ) also found that such superdense 
ETGs are extrem e ly rar e at 2; « 0. On the other hand, 
IValentinuzzi et al] (|2010l ) claim to find a small number of 
compact massive ETGs in the nearby X-ray-selected clusters 
and claim that these galaxy's properties correspond to the 
high-z compact massive ETGs. This last example notwith- 
standing, it is clear observations indicate that most of the 
compact massive ETGs at high-z must increase their sizes 
and must decrease their velocity dispersion from 2: ~ 2 to 
z — without recent star formation, in order to keep their 
old stellar populations. This problem is called as "size evo- 
lution problem of ETGs" . 

Dry (i.e., gas-poor, coUisionless) mergers are a pop- 
ular scenario for the size evolution of the ETGs. Dry 
mergers can increase the galaxy's size more than wet 
(i.e., gas-rich) mergers because of the a bsence of energ; 
dissipation via the gas com ponent (e.g. ICox et al.1 I200' 
ICiotti. Lanzoni. fc Volonterill200"7l ). Furthermore, dry merg- 
ers keep the old stellar populations of progenitor galaxies 
without the addition of new star formation. T here is di- 
rect observational evidence for dry mergers (e.g. iBell et al] 
I2OO6I). Using a semi-analyti c approach to galaxy forma- 
tion, iKhochfar fc Silkl (|2006al ') show that massive ETGs ex- 
perience dry mergers at low redshift {z < 2) and sug- 
gest that dry mergers are the main driver for the size 
evolution of ETGs. There are alternative scenarios to dry 
mergers that involve expansion of the stellar s ystem due 
to a significant mass loss, either by q uasars (|Fan et al.l 
l2008l. iRagone-Fi gueroa fc Granato l201lf ) or stellar winds 
i Damianov et allboogf T However, the amount of the stel- 
lar mass loss is too small to explain the size evolution 
l|Damianov et al]|2009l ') , the compact massive ETGs already 
are quiescent in star formation activity, and they do not 
seem to have sufficie nt gas to puff up the stellar system 
l|Bezanson et al.ll2009l '). 

Dry minor mergers are a more effective process 
for the size evolution of ET Gs than dry major m erg- 
ers. Using the virial theorem , Bezanson et al.l (|2009t ) and 
iNaab. Johansson, fc Ostrikej (|2009l ) make a simple estimate 
of the size growth for dry mergers. They use initial galaxy 
models that consist purely of stellar systems and only as- 
sume energy conservation during mergers on parabolic or- 
bits. They obtain that for dry minor mergers, the galaxy 
sizes increase by the square of the mass increase rather than 
linearly for equal-mass mergers. Semi-analytic models of hi- 
erarchical galaxy formation show that minor mergers are an 
important cha nnel of mass assembly for massive galaxies i n 
z < 1 - 2 (e.g. lKhochfar fc Sil^l2006bl . IGuo fc Whitell2008h . 
Observational studies suggest that even for massive ga laxies 
the major merger rate is low (e.g. iBundv et al.|[2009l '). It is 
not likely that the majo r mergers explain t he extent of the 
size evolution of ETGs (|Tavlor et al.[|2010l '). Therefore, dry 
minor mergers are likely to be the main driver for the size 
evolution of the massive ETGs. 

While a number of numerical simulations have 



been done to s tudy the size evolut io n problem b y dry 
mergers (e.g. iHopkins et al.l |2009|, iNaab et al.l I2OOE , 



Nipoti. Treu. 



Hopkins et al 



BoltonI 
2010d) . 



l2009al . 

they 



iNipoti et ai] l2009b 

insufficient to ex- 
plain the observed difference between high z and local 
ETGs. Several studies have performed high-resolution 
N-body simulations of dry major an d minor mergers to 
investigate the size e volution problem ( Nipoti et al.l 2009al . 



Nipoti et al 



Nipoti et al 



2009b 



, iNipoti et all 12012 . iHilz et al.1 120121 ') 



(2009b|) simulate a set of dry major and minor 



mergers and examine consistency of their results with the 
stellar-mass scaling relations of the local E TGs obtained 
from the Sloan Lens ACS (SLAGS) Survey (jBolton et al.1 
l2008h . They conclude that typical local massive ETGs 
can not have assembled more than ~ 45 % of their stellar 
mass nor grown more than a factor of ~ 1.9 in size via dry 
mergers. While they did explore the effect of minor mergers, 
they limited their study to only the cases of "multiple" and 
"simultaneous" minor mergers. However, theoretical studies 
on merger rates indicat e that simultaneous minor mergers 
arc extremely rare ( e.g. iGuo fc White! l200d IStewart et al.l 
I2OO9, Ho pkins et al.l l2010al ). In the hierarchical structure 
formation scenario in the ACDM cosmology, "sequential" 
minor mergers are more common than simultaneous events. 

In this paper, we investigate the effects of dry major and 
minor mergers on the size evolution of the ETGs with N- 
body simulations, particularly focusing on effects of the se- 
quential dry minor mergers and of the effect of compactness 
of the satellite galaxies on the size evolution. To explore the 
effects of the mergers in inner regions of the compact massive 
ETGs, we take care to simulate dry minor mergers with suf- 
ficient spacial and time resolution. We show that the sequen- 
tial minor mergers of the compact satellite galaxies lead to 
the most efficient size growth and the most efficient decrease 
of velocity dispersion. Furthermore, we construct the merger 
histories of candidates of the high-z compact massive ETGs 
using the Millennium Simulation Database -a public galaxy 
catalog with a semi-ana l ytic model based on the ACDM co s- 
mology (|Springel et al.1 (2005bl . iDe Lucia fc BlaizotI \200t ]- 
and use this to model size growth of galaxies from dry minor 
mergers. Based on the size growth efficiency derived by our 
N-body simulations, we discuss the size growth of galaxies 
from 2 ~ 2 to 2 = 0. 

The paper is organized as follows. In §2, we show the 
importance of minor mergers for evolution of the high-z mas- 
sive (stellar mass : M. ^ IO^Mq) galaxies using the Mil- 
lennium Simulation Database. In §3, we describe the outline 
of the simulations, presenting our initial models and choice 
of simulation parameters in more detail. Results of the anal- 
ysis of the merger remnants of our simulations are presented 
in §4. In §5, we compare our results with observations, and 
discuss the size growth of the compact massive ETGs. In §6, 
we give a summary of the paper. 



2 COSMOLOGICAL MERGER HISTORIES OF 
MASSIVE GALAXIES 

2.1 Sample selection 

To understand the evolution of the ETGs, we examine 
mass assembly histories of the high-z compact massive 
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ETGs. We get merging history information of the dark mat- 
ter haloes and those of the stellar bulges which formed 
in centres of the dark matter haloes, from the results 
of a semi -analytic model in the Millennium Simulatio n 
Datab ase JSpringel et"alll2005bl. JPe Lucia fc BlaizotllioOTi ). 



While iDe Lucia fc Blaizod (|2007l l use strong supernova and 
active galactic nucleus (AGN) feedback models to repro- 
duce observed properties of the brightest cluster galaxies 
(BCGs), they also adjust their model parameters in order to 
be in good agr eement with the obse rvations of the local Uni- 
verse (see also ICroton et alj (|2006h using a similar model) . 
Thus, we use this database as a galaxy formation model 
for all galaxies. We extract galaxy samples at z = whose 
main progenitors have a bulge mass more than IO'^^Mq at 
z — 2.07. These galax ies are candidates for the high-z com- 
pact massive ETGs (jvan Dokkum et"al] l2008h . We found 
1045 galaxies satisfying this condition in the Millennium 
Simulation Database. In the following subsections, we anal- 
yse the merger histories of these sampled galaxies. 



2.2 The important role of minor mergers in 
galaxy mass growth 

Fig. [T] shows the number of the sample galaxies at different 
stellar mass taken at z = (solid lines). These galaxies 
have experienced major and minor mergers from z = 2.07 
to z = 0. We define the minor mergers as mergers with mass 
ratios M2/M1 < 1/4, otherwise we define them as major 
mergers. We also define 'minor merger-dominated galaxies' 
as galaxies for which the minor mergers are the dominant 
form of mass increase during the period z — 2.07 to 2 = 0, 
and we define 'major merger-dominated galaxies' as galaxies 
for which the major mergers are the dominant effect. In the 
left panel of Fig. [1] we show the fraction of the minor merger- 
dominated galaxies in each bulge mass bin, while in the right 
panel of Fig. [T] we show the fraction in each mass bin for the 
main dark matter haloes within which the galaxies reside 
(we call these the "FOF" halo). We find that ~ 65% of 
the total sample galaxies are the minor merger-dominated 
galaxies, and ~ 38% of the total sample galax;ies have never 
experienced major mergers from z = 2.07 to z = 0. 

In the left panel of Fig.[2l we show that fractional cumu- 
lative mass growth of the sample galax;ies by minor mergers 
from z = 2.07 to 2; = 0. Here, the fractional cumulative mass 
growth is the fraction of mass increased by minor merg- 
ers from z = 2.07 to z = to the bulge mass at 2: = 0. 
We find that on average the fractional cumulative mass 
growth by the minor mergers with mass ratios M2/M1 < 1/4 
{M2/M1 < 1/10) is ~ 0.3 (~ 0.2), and more massive galax- 
ies have larger fractional cumulative mass growth by the 
minor mergers. Thus, minor mergers are one of the domi- 
nant processes for the mass growth of our sample galaxies. 
In the right panel of Fig. [21 we show the dependence of 
the fractional cumulative mass growth on the mass of the 
FOF haloes. We find that many galaxies have increased their 
mass significantly by the minor mergers in more massive 
FOF haloes. Therefore, we expect that these minor merg- 
ers increase galaxy sizes, in particular, in high density en- 
vironments. ETGs are likely to be in such high density en- 
vironme nts, which is we ll known as the morphology-density 
relation (|Dressleij[l980l ). For this reason, in this paper we 



concentrate the minor mergers as the one of the important 
process driving the size evolution of ETGs. 



2.3 Stellar merger histories of sample galaxies 

We analyze the merging events of the 1045 sample galax- 
ies, and show that the minor merger events are sequential 
processes rather than simultaneous events between z = 2.07 
and 2 = 0. We obtain merging epochs at which a satellite 
(less massive) galaxy merges into the primary (more mas- 
sive) galaxy from t he results of the semi-analytic model (see 
ICroton et al]l2006l ). We identify 19605 minor mergers, where 
we count multiple minor mergers as one minor merging event 
if they occur in the time interval of ~ 0.35 Gyr, correspond- 
ing to the dynamical time-scale of the dark matter haloes 
of the galaxies at 2 ~ 2. Except for one case, the sum of 
the accreted mass through the minor mergers during this 
time interval is less than 50% of the mass of the main pro- 
genitor galaxy. Therefore, almost all minor mergers occur 
sequentially. 

Fig. [3] shows the stellar merger histories of two sample 
galaxies. One is a central galaxy of a FOF halo which mass 
is 1.5 X 10^^ Mq, and another is a member galaxy of a FOF 
halo which mass is 1.8 x lO^^M©. These two galaxies have 
experienced sequential minor mergers, but have not experi- 
enced major mergers since 2 = 2.07. We find that the total 
merged mass through minor mergers is comparable to the 
stellar mass at z — 2.07. Therefore, the size evolution is 
driven significantly by this process. In this paper, we focus 
such a "sequential" minor merging processes. 



3 N-BODY SIMULATIONS 

In this section, we describe our initial galaxy models, model 
parameters, N-body simulation methods, and analyzing 
methods of merger remnants. 



3.1 Initial galaxy models 

Our model galaxies consist of two components: a stellar 
bulge and a dark matte r halo. We model the dark matter 
mass distribution with a iHernauistI (|l99d ) profile which has 
the total mass, Mdm- 



Pdm{r) 



a.dn 



2ti r{r + Udm)^ ' 



(1) 



scaled to m atch the NFW profile found i n cosmological 
si mulations |Navarro. Frenk. fc W hite 'l997|) , as described 
in lSpringel. Di Matteo. fc Hernauist (,2005a| ). In this model, 
the total mass, Mdm, corresponds to the mass within the 
virial radius of the NFW model and the inner density pro- 
file of the Hernquist profile is also equal to that of the NFW 
model. Satisfying the above condition, the scale radius of 
the Hernquist profile Udm , is related to the scale length of 
the NFW halo and halo concentration Cyir{cvir = r-uir/rs) 
as 



0.dr 



= r-s-\/2[ln(l -f C„ir) — Cyir/{1 + Cvir 



(2) 



We derive the virial radius of the halo from a spherical 
collap se model, r„ir ~ rvir{Mdm, z) (e.g. iBarkana fc Loebl 
l200lh . where 2 is the formation redshift of the halo. Using 
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these parameters of mass, redshift and halo concentration 
for a dark matter halo, we can obtain the corresponding 
halo model for N-body realization. 

We also model the initial distribution of the bulge with 
a Hernquist profile, 

, . h 'L a 

with total stellar mass A/,, quarter-mass radius ri/4 = a 
and half-mass radius ri/2 = (1 + \/2)a- Observationally, a 
galaxy size is usually characterized by its effective radius, 
defined as the radius of the isophote that encloses half of 
the total stellar luminosity in projection. In our simulations, 
we adopt the effective radius Re as the projected radius 
enclosing half of the stellar mass. Here we assume a con- 
stant stellar mass to light ratio. From Eqn. Q, the effec- 
tive radius is related to the quarter- a nd half-mass radii by 
Re = 1.8153ri/4 = 0.752ri/2 (see also iBinnev fc Merrifieldl 
1 19981 ) . We remove particles of the two-component Hernquist 
profile outside r ~ 20adm in our numerical realizations in 
order to avoid the occurrence of a very small number of 
very distant particles. This choice of outer radius does not 
affect the half-mass radius significantly. If the total stellar 
mass and the scale radius are specified, the Hernquist profile 
determines a unique structure. 

We construct the bulge and the dark matter halo to be 
in equilibrium in the total gravitational potential by an N- 
body method. We assume that both components are both 
spherical and isotropic. The particle positions for each com- 
ponent are initialized from the density profile, while the par- 
ticle velocities are drawn based o n each compone nt's equilib- 
rium distribution function fi{E) (jBinnev fc TremainQi200^ . 



HE) 



dV'2 /B^' 



(4) 



where E is the relative energy of each particle of compo- 
nent i, pi is the density profi le of compone nt i, and 1/) is the 
total gravitational potential. ICiottil (| 19961 ) derives the ana- 
lytical distribution function of the two-component Hernquist 
profiles. We use this formula to obtain the particle velocity 
distribution. 

Here, we only briefiy summarize o ur N-body realization 
method, which is described in detail in lKuiiken fc Dubinskil 
(|l994h . A particle's position is first determined by sampling 
from the density distribution. Next the particle's velocity 
vector is randomly sampled from inside a velocity-space 
sphere whose radius is the local escape velocity, Vesc = 
V— 2^. Finally, based on the distribution function of the 
bulge or halo component, we determine whether we select 
the velocity vector or not using the acceptance-rejection 
technique. 

3.2 Model parameters 

For the progenitor galaxies in major and minor merger sim- 
ulations, we assume four galaxy models: a massive ETG, 
two less massive compact ETGs, and a less massive diffuse 
ETG. Parameters of these galaxy models are summarized 
in Table [1] We assume that the massive ETG (model A) 
has the stellar mass M, = lO^^M©, with an effective ra- 
dius Re = l.Okpc, which corresponds to the scale radius 
a = 0.551kpc. These properties are consistent with that 



from observations o f the high-z compact massive ETGs (e.g. 
iTruiillo et al] |2007| ). They are embedded in dark matter 
haloes of Mdm = 10^^ Mq and r^ir at z — 2.0. In this paper, 
we assume dark matter to stellar mass ratio A'Idm/M, — 10 
in all models. We also assume that for all initial dark mat- 
ter haloes, Cvir = 10. Many cosmological simulations pre- 
dict that the halo con centration depends on the mass and 
the formation redshift jBuUock et al.ll200ll '). but we select a 
constant value in favour of simplicity in the model. 

We assume three different satellite galaxies: compact 
(model B), diffuse satellites (model C), and very compact 
satellites (model D). We assume that the formation epoch 
of model B is the same as that of model A, both occurring at 
z = 2.0, and that these models are expected to have the same 
average density. Based on this assumption, we construct the 
stellar system for model B by simply scaling down that of 
model A, given the mass fraction (1/10 of model A) and 
scale radius (1/ \/T0 of model A). We also assume the dark 
matter halo of model B is Mdm = 10^^ M©. For model C 
galaxies, we assume that they have the properties of the 
local ETGs. For their effective radii, we use the observed 
local stellar mass- size relation of early-type galaxies in SDSS 
jShen et al.ll2003h . 



Re = 4.16 



10" M, 



kpc. 



(5) 



The model C galaxies have dark matter haloes with Mdm = 
W^^ Mq and r^ir aX z — 2.0. Recent observations indicate 
that early-type galaxies follow a stellar mass-size relation 
log(i?e) = 5-|-riog(M*) with T ~ 0.6-0.7, not significantly 
dependen t on redshift, and S i ncreasing for decreasing red- 
shift (e.g. iNewman et al.llioi^ '). From this, we assume that 
model A and model D lie on the same stella r mass -size re- 
lation with T = 0.57, as in iNewman et al] (|2012l ). giving 
model D satellites that are more compact than both model 
B and C. We examine effects of the stellar densities of the 
satellites on the properties of merger remnants using this set 
of satellite models. 



3.3 Simulations and test run 



We use the GADGET-2 code (|Springell [ioosi ) . an efiicient 
parallelized treecode, to perform our simulations. We set the 
tolerance parameter a — 0.001 for the relative cell-opening 
criterion in the tree algorithm, the timestep tolerance pa- 
rameter 77 = 0.005, and the force softening e = 0.03 kpc 
for both star and dark matter particles. With this choice, 
the minimum timestep is < 0.02Myr. In this study, we use 
6.0 X 10'' dark matter particles and 6.0 x 10^ stellar particles 
for model A. For model B and model C, we use 6.0 x 10* 
dark matter particles and 6.0 x 10^ stellar particles. Each 
halo and stellar particle has the same mass. 

We should draw attention to the effect of two-body re- 
laxation, since the dynamical time of the compact massive 
ETGs is relatively short compared with the galaxy merg- 
ing time-scale, and we need the simulation time to be much 
larger than the dynamical time of the stellar system. We 
have confirmed that our two-component galaxy model is sta- 
ble at all radii r > a = 0.551kpc for ~ 3.6 Gyr of the simu- 
lation. We define the half-mass dynamical time of a stellar 
system as 
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Table 1. Initial galaxy models for dry merger simulations. First column: name of the initial galaxy model, Model 
A is the compact massive ETG seen at high-z. This is the primary galaxy in our minor merger simulations. 
Model B is the compact, less massive ETG. This is used as a compact satellite. Model C is the less massive 
ETG that obeys the stellar mass-size relation in the local Universe. Model D is the most compact satellite (see 
i|3.2ll . From column 2 to column 10: dark matter halo mass; redshift; halo concentration; virial radius of the 
halo measured at each redshift; number of dark matter particles; stellar mass; effective radius; number of stellar 
particles; softening length. 



Model 




z 


c 


r„ir(kpc) 




m4Mq) 


i?e(kpc) 




e(kpc) 


A 


1012 


2.0 


10 


105 


6.0 X 10^ 


10" 


1.0 


6.0 X 10'' 


0.03 


B 


10" 


2.0 


10 


48.8 


6.0 X lO-i 


lOlo 


0.464 


6.0 X 10^ 


0.03 


C 


10" 


0.0 


10 


120 


6.0 X 10-1 


lOio 


1.15 


6.0 X 10^ 


0.03 


D 


10" 


2.0 


10 


48.8 


6.0 X 10* 


IQio 


0.148 


6.0 X 10^ 


0.03 



Table 2. Physical properties of the merger remnants within the fixed radius, 60kpc, in each run. The properties of the initial galaxy are 
defined within 60kpc. First column: name of the simulation, 2A: major merger of model A, lA: the primary galaxy, 1/5/10: the number of 
minor mergers, B/C/D: type of the satellite galaxy, sq/sm: sequential or simultaneous minor mergers. From column 2 to column 10: stellar 
mass; virial radius of the dark matter halo; half-mass radius of the stellar system; effective radius; average density within Ikpc; average 
density within the effective radius; velocity dispersion; size growth efficiency; efficiency of the change of velocity dispersion (see i|4.5l l. 



Name 


stellar mass 


^vir 




Re 


P(< Ikpc) 


p(< Re) 




a 


7 




(Me) 


(kpc) 


(kpc) 


(kpc) 


(lOiOA/0kpc-3) 


(lO^^Mokpc-^) 


(km/s) 


{Re oc Mf ) 


(ct oc M~'') 


A (initial) 


9.82 X lO*'' 


257 


1.29 


0.968 


0.992 


1.07 


262 






2A 


1.91 X lO" 


309 


3.81 


2.81 


0.933 


0.0864 


253 


1.60 


0.0528 


lAlB 


1.04 X 10" 


261 


1.50 


1.12 


0.961 


0.737 


256 


2.60 


0.435 


lAlC 


1.02 X 10" 


260 


1.40 


1.05 


0.981 


0.878 


258 


2.37 


0.468 


lAlD 


1.05 X 10" 


262 


1.55 


1.16 


0.935 


0.667 


253 


2.69 


0.501 


lA5Bsq 


1.32 X 10" 


275 


2.90 


2.15 


0.852 


0.137 


237 


2.72 


0.350 


lA5Csq 


1.16 X 10" 


268 


1.92 


1.44 


0.951 


0.400 


249 


2.44 


0.310 


lASDsq 


1.35 X 10" 


276 


2.77 


2.08 


0.696 


0.149 


223 


2.38 


0.497 


lAlOBsq 


1.49 X 10" 


301 


4.13 


3.00 


0.853 


0.0581 


226 


2.71 


0.351 


lAlOCsq 


1.19 X 10" 


286 


2.04 


1.54 


0.951 


0.339 


247 


2.42 


0.318 


lAlODsq 


1.70 X 10" 


302 


5.14 


3.76 


0.677 


0.0339 


206 


2.46 


0.440 


lAlOBsm 


1.60 X 10" 


294 


3.78 


2.78 


0.873 


0.0767 


241 


2.17 


0.175 


lAlOCsm 


1.54 X 10" 


287 


3.88 


2.86 


0.964 


0.0706 


243 


2.42 


0.166 



Table 3. Physical properties of the merger remnants within rtrunc = 30, 60, and 90kpc in binary 
minor merger simulations. Each column is similar to Table [2] 



Name 


stellar mass 


Re 


Ue 


a 


7 




(Mo) 


(kpc) 


(km/s) 


{Re oc M^) 


{a oc M-'<) 


fana = 30kpC 


lAlB 


1.01 X 10" 


1.06 


257 


2.71 


0.535 


lAlC 


0.984 X 10" 


0.990 


259 


2.63 


0.758 


lAlD 


1.02 X 10" 


1.11 


254 


2.79 


0.569 


fana = 60kpC 


lAlB 


1.04 X 10" 


1.12 


256 


2.60 


0.435 


lAlC 


1.02 X 10" 


1.05 


258 


2.37 


0.468 


lAlD 


1.05 X 10" 


1.16 


253 


2.69 


0.501 


fana = 90kpC 


lAlB 


1.05 X 10" 


1.14 


256 


2.55 


0.408 


lAlC 


1.03 X 10" 


1.07 


257 


2.28 


0.457 


lAlD 


1.06 X 10" 


1.18 


252 


2.67 


0.534 
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idyn — 



3tv 



16Gp 



1/2 



' 1/2 

2GM* 
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where pi/2 — 3M, /Svrrjyj is the mean density inside the 
half-mass radius ri/2- Here, we consider only the stellar com- 
ponent for simplicity. We estimate the two-body relaxation 
time-scale tr from the following expression, 



tr{r) 



O.IN , 

IniV 



tdyn t 



(7) 



where A*' is the number of stellar particles in a radius r. 
The half-mass dynamical time of model A is tdyn — 5.2 Myr 
and the relaxation time is tr — 2 Gyr which is comparable 
to the total computational time of each of our simulation 
sets. To check whether our two-component galaxy model is 
stable we simulated a single model A galaxy for ~ 3.6 Gyr 
as a test run. Fig. |4] shows several snapshots of the density 
profiles of the stellar system. This figure clearly shows that 
our model is stable in all radii r > 0.2kpc. We checked the 
conservation of the total energy E and the virial ratio {V = 
2T/W), obtaining \AE/E\ < 1 percent and |AV/V| < 0.35 
percent. We also checked the stellar half-mass radius, the 
stellar effective radius and the stellar velocity dispersion at 
3.6Gyr, and we confirmed that changes of both the half-mass 
radius and the effective radius are less than 3 percent, and 
the projected velocity dispersion is nearly constant. These 
changes are smaller than that due to mergers in all runs (see 
Table |2]). In order to avoid the influence of numerical error, 
we stopped our simulations at 3.6 Gyr. 

3.4 Production run set 

We have run N-body simulations of dry major and minor 
mergers to investigate the effects of mergers on the evolu- 
tion of size and velocity dispersion in ETGs. In this paper, 
we mainly focus on the dry minor merger cases. We also 
simulate a dry major merger for comparison. We describe 
the setup of our simulations below. Table [2] summarizes the 
parameters used in our run set. 

Merger orbits are critical to the properties of the 
merger remnants. We assume head-on parabolic mergers 
wi th zero orbital angular momentum in al l our simulations. 
In lBovlan-Kolchin. Ma. fc Quataerd ()2006l ) , who studied dry 
major mergers of early-type galaxies, merger remnant prop- 
erties are highly influenced by orbital parameters of the 
mergers. They showed that the stellar bulges of the merger 
remnants in the cases of radial orbits are larger than in the 
cases of circular orbits. Therefore, our simulations are opti- 
mal cases for the size growth of galaxies. 



3.4-1 Major mergers 

In run 2A, we simulate a major merger for the parabolic and 
head-on encounter of same initial galaxies of model A. The 
initial separation between the two galaxies is 200kpc, and 
the initial relative speed is 308km/s. The merger appears to 
be mostly completed by 2.3 Gyr. 



3.4- S Minor mergers 

We simulate two cases of minor mergers, involving simul- 
taneous minor mergers and sequential ones. We investigate 



how the merger remnants are affected by a different types 
of minor mergers. In each case, we assume three different 
model galaxies (model B, model G, and model D, as in Ta- 
ble [l]). Initially, we assume the central galaxy of model A 
(the "primary galaxy" ) is surrounded by satellite galaxies of 
model type B, C or D. In the simultaneous minor mergers 
case, the satellites are randomly distributed around the pri- 
mary galaxy with the relative distance between each satellite 
and the primary galaxy in the range [130, 170] kpc. During 
the merging events, the satellites are accreted at the same 
time. The merger orbit of each satellite is assumed to be 
parabolic and head-on in the absence of the other satellites. 
We calculate the relative velocity using the total mass of the 
progenitors before removing very distant particles. 

As shown in §2, in hierarchical structure formation sce- 
nario in the AGDM cosmology, there are many sequential 
minor mergers. We study the effects of this sequential minor 
merging in our next case. Here, each satellite merges with 
the primary galaxy from random directions every 0.2 Gyr. 
We simulate five sequential (run IA5Bsq, f A5Csq, IA5Dsq) 
and then ten sequential (run f Af OBsq, lAIOGsq, f Af ODsq) 
minor mergers. In these cases, the period of the sequen- 
tial minor mergers roughly corresponds to the time interval 
from z = 2 to z = 1. We assume the three types of satellites 
(model B, model C, model D), and compare the properties 
of the merger remnants. In the sequential minor mergers, the 
mass ratio of the first merger M2/M1 — 1/10, that of the 
second merger is 1/11, and so on. The mass ratio can reach 
~ 1/20 in the last merger. Thus, the typical mass ratios 
of the minor mergers are 1/20 < M2/M1 ^ 1/10. We also 
simulate single minor mergers (run lAlB, lAlC, lAlD) to 
compare these results with the sequential minor mergers. In 
these test simulations, we assume the initial relative velocity 
between progenitors using the total mass of the progenitors 
after we have removed very distant particles to ensure that 
the satellite orbits are strictly parabolic and head-on. 

Numerical simulations were run on Cray XT4 at Genter 
for Gomputational Astrophysics, CfCA, of National Astro- 
nomical Observatory of Japan. 



3.5 Definition of merger remnants 

In this subsection, we describe our analyzing methods for 
the merger remnants. In the major merger simulation (run 
2A), the final state is defined when the global virial ratio of 
the system reaches an almost constant value. In the case of 
the minor mergers, it takes more time for the minor mergers 
to attain complete relaxation, since the dynamical friction 
time-scale of the minor mergers are longer than the merging 
time-scale of the major merger. Therefore, as described in 
H3.3I we stop our runs a.t t = 3.6 Gyr in the minor merger 
simulations, and analyze the merger remnants. At this time, 
there are several surviving cores of the satellite stellar sys- 
tems in the remnants from the 10 minor mergers. However, 
the mass of each core is less than 5 per cent of the mass 
of the merger remnant which we define below. We exam- 
ine the effects of these surviving cores on properties of the 
merger remnants in M.31 and find that these cores have little 
effect on the size growth efficiency (defined in Sj43]) of the 
remnants. 

We analyse the remnant particles at the end of the sim- 
ulations. For the dark matter particles, we define the merger 
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remnants as the particles which are gravitationally bound. 
For the stellar particles, it is difficult to use that criterion, 
since in the minor merger simulations, the central galax- 
ies are surrounded by very diffuse stellar particles beyond 
> 50 kpc that are gravitationall y bound by the dark mat - 
ter haloes of the remnants. Since Ivan Dokkum etlo] (|201(]| ') 
use 75kpc as the robust galaxy radius to define total stellar 
mass in their stacked images, and determine the effective 
radii, we use a radius rtrunc = 60kpc to define the stellar 
remnants, with our reasoning described below. We call the 
stellar mass within rtrunc the remnant mass. In r > rtrunc, 
the stellar surface density of the remnants is too low to 
be observed. In fact, the surface brightness difference be- 
tween O.lkpc and 60kpc are A/i > 13mag arcsec"^ for the 
remnants, assuming a constant mass-to-light ratio with ra- 
dius. Even state-of-the-art observations cannot observe such 
a large range of mag nitudes (Truiillo ct al. 20li). We exam- 
ine various radii, rtrunc, to find a suitable rtrunc for defin- 
ing the remnants. We have checked for ''"trunc = 30, 
60, and 90kpc in the binary minor merger simulations as 
shown in Table [3l We find that Re slightly increases from 
rtrunc = SOkpc to rtrunc = 60kpc, and is almost constant be- 
tween rtrunc = 60kpc and rtrunc ~ 90kpc. Thus, we choose 
rtrunc = 60kpc in the following analysis of merger remnants. 

We analyse the three-dimensional and projected prop- 
erties of the stellar merger remnants. The stellar half-mass 
radius is the radius of the sphere within which half of 
the mass of the remnant stellar system is enclosed. We 
adopt the effective radius Tie as the projected radius en- 
closing half of the remnant mass. We also define line- 
of sight stellar velocity dispersion of the remnants as in 
iBovlan-Kolchin. Ma. fc Quataerd ()2005l l . We use the surface 
mass density weighted velocity dispersion within R^, 



j^; aLiRMR)RdR 

J^^^ T,{R)RdR 



(8) 



where is the surface mass density and aios is the line- 

of-sight velocity dispersion. We use 2e as the lower limit of 
the integrations in Eqn. ([8]), where e is the softening length. 
Taking into account the projection effects, we average Re 
and (Te in 100 random projections. 

To derive major and minor axes of the merger remnants, 
we obtain the principal axes of the merger remnants using 
the reduced moment of inertia te nsor (jPubinski fc Carlbergj 
ll99ll . lBovlan-Kolchin erahlbOOSD . 



mining these principle axes has to be done iteratively, start- 
ing with the assumption of spherical symmetry, q = s = 1. 



4 RESULTS 

We have run N-body simulations of the dry major, the si- 
multaneous minor and the sequential minor mergers as de- 
scribed in ij3l Physical quantities of the merger remnants are 
summarized in Table [2l 



4.1 Snapshots 

Figure [S] shows five snapshots for three runs, lAlOBsm, 
lAlOBsq and lAlOCsq (Table [2]). Elapsed time from the 
beginning of each run is written in the box of each snap- 
shot. As described in Sj3l we stop runs after t = 3.6 Gyr 
of evolution in the minor merger simulations, and analyze 
the merger remnants nt t — 3.5 — 3.6Gyr, to allow a small 
number of surviving satellites to be among the remnants. 

In this final state, there are the several surviving cores 
of the satellite stellar systems among the remnants of the 10 
minor mergers. However, most of stars and the dark mat- 
ter component of the satellite galaxies have already been 
stripped or merged into the primary galaxies. We show that 
these surviving cores have only a small effect on properties 
of the merger remnants in i ]4.3l In the major merger (run 
2A), the first crossing of galaxies is at ~ 0.5Gyr, and the 
merger appears to be completed by 2.3Gyr. 



4.2 Density profile 

One of the most revealing properties of the remnants is their 
density profiles. Fig.[6]shows the angle-averaged stellar den- 
sity profiles of the merger remnants for five runs. In this 
figure, we show three types of final density profiles for each 
run: the profile of the stars that initially belong to the pri- 
mary galaxy (we call "primary stars"), that of the stars that 
initially belong to the satellite galaxies (we call "satellite 
stars"), and that of all the stars. For comparison, the den- 
sity profile of the initial galaxy model is also shown. The 
density profiles for all the stars differs significantly between 
each run. Inside the inner region of the half mass radii (sum- 
marized in Table [2]), the gravity of the stellar mass is closely 
related with the velocity dispersion of the stars in this re- 
gion, as shown in i]4.4l Below, we discuss these differences 
in each of the inner and outer sections of the remnants. 



N 



(9) 



(10) 



where ai is the elliptical radius and q and s are the axis 
ratios with s ^ g ^ 1. q and s are calculated through 



/ -^yy 



r xi/2 /r xi/2 

and s = i — 



(11) 



where Ixx, lyy and Izz are the principal component of the 
tensor. A'' is the number of particles in the remnants defined 
by all stellar particles within rtmnc- The procedure in deter- 



4-2.1 Inner region of the stellar system 

To make clear the difference between the five runs in the 
inner regions, we show the cumulative mass distributions of 
the star particles of the remnants in Fig.[71 We first focus the 
distributions of the primary stars. Comparing the distribu- 
tion of the primary stars of the initial state with that of the 
final state for the minor mergers of the compact satellites 
(run lAlOBsm, lAlOBsq, and lAlODsq), we show that the 
decrease in the density profiles is due to expansion of the pri- 
mary stars. Such expansion is caused by dynamical friction 
heating from the satellites (lEl-Zant. Shlosman. fc HoffmanI 
l200ll , lMa fc Bovlan-Kolchinll2004 ). We see in Figure [71 that 
this expansion is smaller for the minor mergers with the 
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diffuse satellites (run lAlOCsm and lAlOCsq). As shown 
in Fig. [5] the compact satellites can survive for a longer 
time-scale than the diffu se satellites in the tidal field of the 
primary galaxy (see also iBovlan-Kolchin fc Mall2007l ). This 
increased lifetime means that the dynamical friction heating 
can be more effective in expanding the primary stars. There- 
fore,the compact satellites are more efficient at heating than 
the diffuse satellites. 

In the inner region, the density profiles of the satellite 
stars in the remnants depend on the compactness of the 
satellites. From Fig. [6] for the compact satellites, we can see 
that there are more satellite stars in the half mass radii than 
for the diffuse satellites. Once again, this is due to the com- 
pact satellites surviving against the tidal field of the primary 
galaxy to reach its inner regions. For the minor mergers with 
the diffuse satellites, the satellite stars have a smaller con- 
tribution to the density profiles in the half mass radii. This 
is because the diffuse satellites are easily disrupted by tidal 
stripping and tidal shocking from the primary galaxy. This 
result is consistent with the minor merger simulation s of a 
single satellite galaxy bv iBovlan-Kolchin fc Mai (|2007l 'l. 

Compared with the sequential minor mergers, simul- 
taneous minor mergers find a higher number of satellite 
stars in the half mass radii. It is known that in a time- 
varying self-gravitational potential, particles exchange their 
energies with each other, altering the energy distribution of 
the particles. This process is known as violent relaxation 
l|Lvnden-Bel]||l967l ') and produces tightly bound particles. 
i]4.4l will show in more detail that the satellite stars of the 
simultaneous minor mergers have experienced this process. 
We will also show in the case of sequential minor mergers, 
that the violent relaxation process is less effective, and the 
number of tightly bound particles is smaller than for the 
simultaneous minor mergers in i ]4.4l 

4-2.2 Outer region of the stellar system 

Fig. |6] shows that the increase in density in the outer re- 
gions of the density profile is dominated by the accretion of 
the satellite stars onto the outside of the half mass radii. 
The density profiles of the primary stars is hardly changed, 
as shown in Fig. |6] and follows a r^^ profile in the outer 
region. On the other hand, the slope of the density pro- 
files of the satellite stars follows an r~^ up to several 100 
kpc. Then for radii larger than the several 100 kpc, the 
slope declines as p oc r"**. These difference show that the 
two components have distinct origins; an already formed 
system and a newly accreted system. In particular, in the 
case of the sequential minor mergers of the diffuse satel- 
lites (run lAlOCsq), the satellite stars contribute to the en- 
velope of the density profile of the remnant as an excess 
co mponent, as shown in Fig. [6] Similar results arc obtained 
bv lAbadi. Navarro. fcSt cinmctz (200&).,Abadi ct al.. (2006 ) 
simulate an isolated disc galaxy with cosmological initial 
conditions and show that the density profile of the accreted 
stellar component follows the slope p oc r~^ near the edge 
of the primary stellar component of the galaxy. In addi- 
tion, the slope reaches p oc or steeper near the virial 
radius of the system. Several recent observations of com- 
pact massive ETGs at low-z (z < 0.5) show that the surface 
brightness profiles of the ETGs are well fitted with two- 
component Sersic profiles jStockton. Shih. fc Larsonir2010l . 



IShih fc Stocktonll201ll ). While these Sersic indexes are rel- 
atively small (n< 2), these galaxies may have experienced 
the size growth by minor mergers like our simulations. 

4-2.3 Dark matter halo 

Fig. \7\ also shows cumulative mass distributions of the dark 
matter particles of the merger remnants of the five runs. 
The mass of dark matter haloes are smaller than the stellar 
mass in the inner regions. However, change of dark matter 
mass distribution is important, since this affects the velocity 
dispersion of the stars in this region, as shown in H4.4I In 
the case of the minor mergers of the compact satellites, the 
dark matter particles that initially belong to the primary 
galaxy (we call "primary DMs") expand due to dynami- 
cal friction heating. Compared with the simultaneous minor 
mergers, sequential minor mergers have only small number 
of dark matter particles that initially belong to the satellite 
galaxies (we call "satellite DMs") contribute to the density 
profiles within lOkpc from the centre of the merger remnant. 
For simultaneous minor mergers, a significant amount of the 
satellites' DM is deposited in the central lOkpc. Thus, se- 
quential minor mergers of compact satellites are the most 
efficient in expanding of the dark matter distribution and in 
decreasing of the central dark matter density. 

4.3 Stellar surface density profile 

We show the angle-averaged stellar surface density profiles of 
the merger remnants in Fig. |8l For the evolution in the sur- 
face density profiles, the choice of projection direction does 
not significantly alter the result. For the major merger (run 
2A), the surface densities in r < Ikpc increases from the 
initial state. On the other hand, for the minor mergers, the 
surface densities are approximately stationary in r < Ikpc. 
In > 2kpc, the growth in the surface densities is dominated 
by the accretion of the satellites onto the outer envelope. 
The latter c hange agrees well wi t h the recent observational 
results from Ivan Dokkum et"al] (120101 1 and is an example 
of "inside-out growth". These profiles depend on the com- 
pactness of the satellite galaxies. In Fig. |8l we show that 
inside-out growth is significant when the satellite galaxies 
are diffuse. 

In the case of 10 sequential minor mergers, there are 
several surviving cores of the satellite stellar systems within 
''trunc of the remnants. Due to these cores, the remnant sizes 
might be artificially large in our analysis. We examine effects 
of these cores in determining sizes and mass of the merger 
remnants by removing the satellite cores from the surface 
density profile of run lAlOBsq. We find that the remnant 
after the removing is ~ 0.67 x 10^" Af© less massive and 
~ 0.32kpc smaller than with the cores included. This size 
change is much smaller than the size growth of the remnant 
of lAlOBsq that increased from Ikpc to 3kpc. Therefore, in 
our analysis, these cores do not affect the size growth of the 
remnants significantly. 

4.3.1 Sersic profile fitting 

In order to compare the structures of our merger remnants 
with observational results quantitatively, we fit the stellar 
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Table 4. The Sersic index. The b — c plane means the pro- 
jection along the major axis. The remnant is projected in the 
plane with the intermediate axis b and the minor axis c. The 
a — b plane is the projection along the minor axis. This plane 
contains the major axis a and the intermediate axis b. 



Name 


b ~ c plane 


a — b plane 


A (initial) 


4.07 


4.04 


2A 


5.00 


6.16 


lAlOBsm 


6.51 


6.32 


lAlOBsq 


11.3 


8.47 


lAlOCsm 


8.41 


8.27 


lAlOCsq 


4.83 


4.76 


lAlODsq 


4.66 


7.87 



surface density profiles of the remnants along the maj or and 
the m inor axes derived by the method in §3.5 with the lSersi"3 
(| 19681 ^ i?'/™ law: 



I{R) = lo exp 



(12) 



where fe(m) ~ 2m - 1/3 + 4/(405m) (jCiotti fc Bertinlll999h . 
We choose the radial range [0.1, 10] kpc for this fitting. In 
this range, the surface density profiles are well fitted with a 
single Sersic profile. We note that the Sersic index m is very 
sensi tive to the fitting range (see also iBovlan-Kolchin et al] 
I2OO5I ). 

We find that the Sersic indexes of the remnants in- 
crease from our initial model in all runs. These val- 
ues are consiste nt with the o b serve d massive ETGs (e.g. 

Del Principi |2002| . iGraham fc GuzmanI 
lAlOBsq, the Sersic index 
of the marked inside-out 



Bertin. Ciotti. 



2OO3I ). In particular, for run 
increases significantly because 



growth. On the other hand, for run lAlOCsq, the Sersic in- 
dex increases less. As described above, this is because that 
the stars of diffuse satellites are easily stripped by the tidal 
field of the primary galaxies, and few satellite stars are added 
inside of lOkpc from the centre of the remnant. The increase 
of the Sersic index significantly depends on the compactness 
of the satellites and the manner of minor mergers (sequential 
or simultaneous). 



4.4 Velocity dispersion profile 

In Figure [9] we show the evolution of the line-of-sight stel- 
lar velocity dispersion profiles for the four runs with minor 
mergers. As with the surface density profiles, the choice of 
projection does not affect our main results. Fig.[9]shows that 
in the cases of minor mergers with the compact satellites, the 
velocity dispersion of the stellar systems in the central region 
decreases from the initial model more strongly than for the 
minor mergers with the diffuse satellites. This is because the 
central region is gravitationally heated by the compact satel- 
lites. We compare the velocity dispersion of the remnant of 
run lAlOBsm with that of run lAlOBsq at ~ 2.4kpc, which 
is near the effective radii of these remnants. It is obvious 
from Figure [9] that for run lAlOBsq, the velocity dispersion 
in this radius is smaller than for run lAlOBsm. This is be- 
cause for run lAlOBsq, few satellite stars and satellite DMs 
are deposited within this radius. For run lAlOBsm, rela- 



tively large amounts of the satellite stars and satellite DMs 
are deposited within this radius, and the velocity dispersion 
of the remnant increases from the initial value in this region. 
Therefore, in the case of the sequential minor mergers of the 
compact satellites, the velocity dispersion is smaller in this 
region. 

To interpret the difference in the mass deposited dur- 
ing both runs, we show the differential energy distributions 
N{E) for the initial and final stellar components in Fig. 
IIOI where N{E)dE represents the number of particles of 
specific energy between E and E + dE. For the simultane- 
ous minor merger case, the distribution of the satellite stars 
broadens during the mergers and the low energy stars are 
more bound than the initial state. This broadening of N{E) 
is simi lar to the behavior of N{ E) of cold collapse simula- 
tions (|Binnev fc Tremaind [ioosi ') . This means that in this 
case the violent relaxation is effective for the satellite stars. 
Due to such evolution, the satellite stars are more bound 
and strongly invade the centre of the primary galaxy. On 
the other hand, there is smaller change of N{E) for the se- 
quential minor merger case. This means that in this cases 
the violent relaxation is not effective and the satellite stars 
can not invade the inner region. In both simulations shown 
in Fig. 1101 the number of the most bound primary stars are 
reduced. This unbinding is induced by the dynamical friction 
heating during the minor mergers. 

We also investigate anisotropy of the velocity dis- 
persions of the stellar merger remnants. We measure an 
anisotropy parameter for the remnants 



(13) 



where ae and ar are the azimuthal and radial velocity disper- 
sions with respect to the centre of the remnants, respectively. 
Here, we assume ag = a^. In the bottom panels of Fig.[9l we 
show the anisotropy of the merger remnants. For all runs, ra- 
dial anisotropy increases with radii. It is obvious that in the 
case of sequential minor mergers, radial anisotropy is smaller 
than that for simultaneous minor mergers. This result is 
consistent with Atlu uiassoula (200i) who shows that radial 
anisotropy becomes smaller for larger time intervals between 
mergers in sequential galajcy mergers in galaxy groups. We 
note that for the sequential minor mergers satellite orbits 
are radial with respect to the 'initial' position of the pri- 
mary galaxy, and actually these orbits may be more circular 
than for the simultaneous minor mergers. This may affect 
the small radial anisotropy for the sequential minor merg- 
ers. Our results suggest that the manner of minor mergers 
influences the radial anisotropy of the merger remnants. 

4.5 Evolution of stellar mass, Re, and ae 

We summarize the results for stellar masses, the effective 
radii Re, and the velocity dispersions ae of the merger rem- 
nants in Table [21 The minor mergers of the compact satel- 
lites are more efficient in increasing the mass of the rem- 
nants than those of the diffuse satellites. More than half of 
the satellite stars in the compact satellites contribute to the 
mass of the remnants. On the other hand, in the case of the 
diffuse satellites (model C), the satellites are disrupted by 
tidal stripping and shocking during the merging, and most 
of satellite stars of the diffuse satellites are distributed in 
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the envelope of the remnants. These are too diffuse to ob- 
served and do not contribute to the mass of the remnants 
(run lAlOCsq). 

Changes of and are influenced both by the manner 
of the minor mergers and by the compactness of the satel- 
lite galaxies. For the minor mergers of the compact satel- 
lites, the sequential cases cause larger size growth and a 
greater decrease in the velocity dispersion compared to that 
in the simultaneous cases. On the other hand, for the minor 
mergers with the diffuse satellites, the simultaneous cases 
cause larger size growth, while the decrease in the velocity 
dispersion is weak in both the sequential and the simultane- 
ous cases. There are two processes which are closely related 
with the changes of Re and a^- The first process is the ex- 
pansion of the primary stars by dynamical friction heating. 
This increases the size of the primary galaxy and decreases 
the velocity dispersion, if it is principally determined by the 
mass distribution of the primary stars. In the case of minor 
mergers with the compact satellites, this expansion is more 
effective in both the sequential and the simultaneous cases 
than for the minor mergers with the diffuse satellites (Fig. 
[7|. The second process is mass deposit of the satellite stars 
by the minor mergers. The mass deposit increases the size of 
the remnant efficiently for all runs. It is important for change 
to the velocity dispersion whether the satellite stars are de- 
posited within the effective radii of the merger remnants or 
beyond it. If a significant number of the satellite stars are 
deposited within the effective radii, the velocity dispersions 
of the remnants increase. In the case of the sequential minor 
mergers, the mass deposit within the effective radii is small. 
In this case, the velocity dispersion is mainly determined 
by the expansion of the primary stars, and the velocity dis- 
persion decreases. These two processes also cause a changes 
in the density and the velocity dispersion profiles (see i\4.2l 
i]4.4[) . We note that comparison of run lAlOBsq, lAlOCsq, 
and lAlODsq shows that for the minor mergers of the com- 
pact satellites, the effective radii increases more than for the 
minor mergers of the diffuse satellites. This result is appar- 
ently contrary to expectations bas ed on the virial t heorem 
and the energy conservation fe.g. iNaab et al.ll2009l ). How- 
ever, because of tidal disruption of the diffuse satellites as 
described above, the size and mass increase is weak in the 
case of the diffuse satellites. 

We present two important relations from our results. 
The first is the relation between the size growth and the mass 
growth. The second is the relation between the decrease in 
velocity dispersion and the mass growth. We parameterise 
the relation between the size growth and the mass growth 



ReJ ^ (MfV 
Re,i \Mi) 



(14) 



to quantify the efficiency of the size growth (hereafter we call 
a the size growth efficiency), where Re,i and R^j are a initial 
and a final effective radius, Mi and Mf are a initial and a 
remnant stellar mass, respectively. We also parameterise the 
relation between the change of the velocity dispersion and 
the stellar mass growth as 



CJe,i \ Mi ) 



(15) 



persion. Table [2] shows a and 7 for each run. As shown in 
Table [2] the cases of the sequential minor mergers, a and 
7 are generally higher than that for the simultaneous minor 
mergers. In particular, we show that the sequential minor 
mergers of model B are the most efficient in the size growth, 
a ~ 2.7. 7 is larger for runs with model D, the most com- 
pact satellite. In the case of the sequential minor mergers of 
the diffuse satellites (run lAlOCsq), a and 7 are also larger. 
In this case, only the expansion is effective. However, here 
the satellite stars do not contribute to the mass increase of 
the merger remnant because of the tidal disruption of the 
diffuse satellites as described above. Therefore, the minor 
mergers of the diffuse satellites do not contribute the size 
growth significantly. 

As described in the previous section, in the case of 10 
sequential minor mergers, there are several surviving cores 
from the satellite stellar systems in the remnants within 
ftrunc. Using the prescription for removing satellite cores 
for lAlOBsq, we find that the size growth efficiency a ~ 2.7 
with or without the satellites core. Therefore in our analy- 
sis, these cores do not affect the size growth efficiency of the 
remnants significantly. 



5 DISCUSSION 

5.1 Size growth efficiency 

As shown in Table [2] the size growth efficie ncies a are 
highe r than for an analytical estimation (e.g. iNaab et al.l 
2OO9I) or found numeric ally in previous works (e.g. 
Bovlan-Kolchin et al.ll2006l l. We check these high efficien- 
cies using a simple analytical argument: Under the assump- 
tions of energy conservation, single component galaxie s, and 
a parabolic orbit in a minor merger, iNaab et all (|2009D show 
that the ratio of the final to initial gravitational radius is 



(1 + n^) 



(16) 



Here, we define 77 = M2/M1 and e = {vf)/{v2), where Mi 
and (wj ) are the mass and the velocity dispersion of the more 
massive galaxy before the minor merger M2 and (wl) are 
those of the less massive galaxy. This formula is also valid for 
two-component galaxy models if the dar k matter h alo profile 
strictly follows the stellar density profile (|Nipoti et al. 201 J). 
We adopt this assumption in the following discussion, rg in 



Eqn. (|16|) can be replaced with Re, where rg and Re satisfy 
the following equation. 



E = -f 



Re 



IGMl 

'2 rg ' 



(17) 



where E is the energy of the stellar system and / is a struc- 
tural parameter that depends on the stellar and dark matter 
density profiles. We obtain 



ReJ _ ff (I+Vf 



Rea fi (1 + r?e) 



(18) 



to quantify the efficiency in the change in the velocity dis- 



where /i and // are the structural parameters of the pro- 
genitor and remnant stellar systems. Here, we assume that 
two progenitors have the same structural parameter /i. This 
assumption is valid for runs lAlB and lAlC. Thus, if / in- 
creases in a dry minor merger. Re increases more than for 
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Table 5. Physical properties of progenitors and merger rem- 
nants in binary minor merger simulations. Note that we define 
the merger remnants as bound particles here. Each column is 
similar to Table [2] except for fifth column, structure parameter 
of the galaxy. 



Name 


stellar mass 


Re 


a 


/ 






(kpc) 






modelA 


9.97 X IQio 


0.995 




0.173 


modelB 


0.996 X IQi" 


0.453 




0.177 


modelC 


0.994 X 10^'' 


1.11 




0.173 


modelD 


0.999 X 10"' 


0.261 




0.163 


lAlB 


10.8 X lO^'' 


1.21 


2.40 


0.185 


lAlC 


10.7 X lO^" 


1.16 


2.20 


0.183 


lAlD 


10.9 X IQi" 


1.23 


2.54 


0.183 



the case where // = fi. From Eqn. (|18p . the size growth 
efficiency can be written as 

rflnM. ln(l + r?) (! + ??)' 

We estimate a from ff/fi, r], and e from our numerical 
results with Eqn. (|19p . In the case of run lAlB a = 2.2, 
since ////i = 1.07, 77 = 0.1, and e = 0.446. For run lAlC 
a = 2.3, since fj/fi = 1.06, r] = 0.1, and e = 0.297. For run 
lAlD a = 2.0, since ////i = 1.06, 77 = 0.1, and e = 0.580. 
We analyse the merger remnants of run lAlB, lAlC, and 
lAlD (see Table [5]), defining the merger remnants as bound 
particles that are consistent with the energy conservation 
argument here. We obtain size growth efficiencies of a = 2.4 
(run lAlB), and a = 2.2 (run lAlC). These values agree 
well with the analytical argument. On the other hand, in 
the case of run lAlD, a = 2.5 which does not agree well 
with the analytical result. In this case, the change in the 
stellar and dark matter halo profiles is larger due to the 
strong dynamical friction heating. Thus, this disagreement 
may stem from the assumption in the analytical argument 
that the dark matter halo profile follows the stellar density 
profile, which is failing to hold in this case. Alternatively, the 
disagreement may be because the assumption that the two 
progenitors have the same structural parameter /i is invalid 
in this case. Therefore, we suggest that our high size growth 
efficiency can be interpreted as the increase of structural 
p arameters / of merger remnan ts (see also a similar analysis 
in lBovlan-Kolchin et al.l (|2006l )). Similar results of the high 
size grow th efficiency for min or mergers on radial orbits are 
shown in lNipoti etM (|2012l ). 

5.2 Comparison with observations 

Our main focus is whether the dry minor mergers can ex- 
plain the size evolution problem of ETGs and reproduce 
their local scaling relations, in particular, the stellar mass- 
size and the stellar mass-velocity dispersion relation. We 
compare our numerical results with these results here. 

5.2.1 Stellar mass-size relation 

In Fig. 111! we compare our numerical results with the lo- 
cal st ellar mass-size relation given by the SDSS (|Shen et al.l 
I2OO3I ). We show that the sequential minor mergers of the 



compact satellites cause the most efficient size growth, while 
all results move the compact massive ETGs toward the local 
relation in Fig. 1111 Interestingly, the sequential minor merg- 
ers have a ~ 2.7 (Table[2]), which is much larger value than 
that obtained vi a the simple analysis based on the viria l 
relation, a ~ 2 (|Bezanson et al.l |2009| . iNaab et al.l 120091 ') . 
with the assumptions of energy conservation, single com- 
ponent galaxies, and parabolic orbit of minor mergers. The 
sequential minor mergers of compact satellites can evolve a 
compact massive ETGs to that with the the average size of 
ETGs in the local Universe by a factor of ~ 2 mass growth. 
As shown in §2, such mass growth is likely to occur for the 
majority of the high-z massive galaxies. 

Our high growth efficie ncy is consisten t with the 
observational constra int by iBezanson et all Hooi). In 
IBezanson et al.l (|2009l ') . they suggest that the size growth ef- 
ficiency, a, must be a > 2 in order to satisfy that the number 
density of descendants of the high-z compact massive ETGs 
does not exceed constraints imposed by the z = galaxy 
mass function. Our result for the size growth efficiency for 
the sequential minor mergers, a ~ 2.7, is consistent with 
thei r results. 

iNewman et"aLl (|2012h estimate the size evolution of 
ETGs using the CANDELS survey and conclude that ob- 
served rapid size growth from 2; ~ 2 to z ~ 1 cannot be 
explained by the dry minor merger scenario, if they as- 
sume Q ^ 1.3 — 1.6 motivated by the previous studies (e.g. 
iNipoti et al.l [2009bl ). However, our results show a signifi- 
cantly higher size growth efficiency a ~ 2.7 which may in- 
validate their conclusion. 



5.2.2 Stellar density 

As shown in Table [21 we find that the mean stellar densi- 
ties within Ikpc, p(< Ikpc), decrease by only a small factor 
in all runs. On the other hand, for the sequential minor 
mergers of the compact satellites, the mean stellar density 
within the effective radius, p(< Re), decreases by more than 
an order of magnitude. This evolut ion is consistent with 
several recent observatio nal studies (jBezanson et al ] |2009l . 
Ivan Dokkum et"al] I2OI0I ) which show evidence to suggest 
that the compact massive ETGs at high-z are the central 
stellar systems in the normal nearby ETGs. This scenario is 
called the "inside-out" growth of ETGs. We note that the 
decrease of p(< Re) in our results exceeds that from simple 
analysis (jNaab et aLll2009l l. corresponding to a size growth 
larger than that of the simple estimate. 



5.2.3 Stellar mass-velocity dispersion relation 

In Fig. ll2l we compare our results with the local stel l ar mass - 
velocity dispersion relation derived bv lNipoti et al.l (|2009bl l. 
We find that the remnant of the sequential minor mergers 
of the compact satellites lies on the local relation, and the 
velocity dispersions of all minor mergers decrease toward 
the local relation. We note that the decrease in the veloc- 
ity dispersion for almost all runs is less effective than that 
derived by simple analysis using the virial relation, 7 ~ 1/2 
(|Naab et al.ll2009l ). 
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5.3 Size growth expected from the Millennium 
Simulation Database 

As shown in §2, we analyze the Mil l enniu m Simulation 
Database given bv lDe Lucia fc BlaizotI (|2007l ') to derive the 
cumulative mass growth via minor mergers. In this subsec- 
tion, we present the size growth for the sample galaxies de- 
rived by the Millennium Simulation Database, using the size 
growth formula given by Eqn. (|14|) assuming the size growth 
efficiency a. We confirm that for ~ 95% (~ 91%) of our sam- 
ple galaxies, ~ 80% (~ 90%) of the bulge mass at z = con- 
sist of stars that belonged to their progenitor galaxies. Thus, 
they mainly inc rease their masses b y dry mergers. This is 
mainly because iDe Lucia fc BlaizotI (|2007l ) use strong su- 
pernova and AGN feedback models. A small fraction of our 
sample galaxies increases their masses by star formation. 
These 'wet mergers' may influence their size increase. In 
this discussion, we neglect these 'wet mergers', since their 
number is small. Here, we divide the sample galaxies into 
two categories: brightest cluster galaxies (BCGs) and nor- 
mal ETGs. We define the BCGs as the galaxies that are 
the central galaxies in FOF haloes with mass of > 10^'^ Mq, 
otherwise we define them as a normal ETGs. 

For the sample galaxies, we derive the average size 
growth factor {Rej/Re,i) by the following procedure: First, 
for each sample galaxy we calculate the mass increase 
through minor mergers from z — 2.07 to z = 0, Amminor. 
Second, we calculate the mass growth factor Mf /Mi by the 
minor mergers from z = 2.07 to « = 0, where Mi is the 
bulge mass at 2: = 2.07 and Mf = Mi -t- Amminor. Third, 
we calculate the size growth factor Re.f/Re.i with Eqn. (|14p 
for each sample galaxy, where we assume a size growth ef- 
ficiency a. Here, we adopt the same size growth efficiency 
for all sample galaxies. Finally, we average the size growth 
factors of all sample galaxies as {Rej /Re,i) ■ 

We show the dependence of the average size growth fac- 
tor on the adopted size growth efficiencies in Fig. 1131 It is ob- 
vious that a > 2.3 is needed to explain the observed size evo- 
lution of ETGs from z ~ 2 to z = (a, factor of three to five 



in size growth, e.g. 


Truiillo et al. 2006. 


Truiillo et al.l 20071. 


Buitrago et al.ll200S 


, Cimatti et al.ll2008 


.Ivan Dokkum et al.l 


2008t). This constraint is consistent with the high size growth 



efficiency derived by our numerical simulations. On the other 
hand, the lower size growth efficiency a = 1.3 d erived by the 
simulta neous dry minor merger simulations of iNipoti et al.l 
l|2009bl ) is insufficient to explain the observational results. 
Therefore, we suggest that in the dry minor merger scenario, 
the sequential minor mergers we performed in this paper are 
responsible for the size evolution of the ETGs. 

We should note that we use the maximum size growth 
efficiency in the above discussion, whereas the choice of 
a ~ 2.7 corresponds to the maximum predicted size increase 
and the assumption that all minor mergers are those with 
mass ratios of M2/M1 < 1/10 and their orbits are parabolic 
and head-on. In fact, there are a wide range of mass ratios 
for min or mergers i n cosm ological merger histories. For ex- 
ample, lOser et "al] (|2012l ) find that the typical mass ratio 
of stellar mergers is 1:5 with cosmological hydrodynamical 
simulations. Additionally, there are satellite galaxies which 
have different properties, e.g. morphology, compactness, and 
gas content, which may infiuence the size growth of com- 
pact massive ETGs. In the future studies, we will take into 



account the mass ratio, orbits, and the many satellite prop- 
erties for more robust predictions based on the dry minor 
merger scenario. 

In Fig. 1141 we show the individual size growth factor 
for all sample galaxies as a function of the bulge mass at 
z = 0. Fig. [15] is the same as Fig. [14] but as a function of 
the FOF halo mass to which the each sample galaxy belongs 
at 2; = 0. In this figure, we assume a — 2.7. We note that 
there are a considerable number of the galaxies which have 
a small size growth factor from 2 = 2.07 to 2: = from 
dry minor mergers (Fig. I14(). These galaxies reside in dark 
matter haloes which have a wide range of mass from galaxy- 
size haloes to cluster-size ones (Fig. I15[l . and almost all of 
them are normal ETGs. In addition, they also do not in- 
crease their size by dry major mergers. In fact, in the galax- 
ies which have a size growth factor less than two, ~ 80% 
are minor merger-dominated galaxies. Thus, these galaxies 
are expected to keep their size from 2 ~ 2. Such galaxies in 
cluster-size haloe s may correspond t o near by compact mas- 
sive ETGs which IValentinuzzi et al.) (|2010l ) claim to find in 
nearby X-ray-selected clusters. 



5.4 Comparison with previous work 

Our simulation results show a high size growth efficiency 
and decrease in the velocity dispersion even in our low- 
est growth efficiency case with simultaneous minor merg- 
ers, compared with results from iNipoti et al] (|2009bl ) . This 
disagreement may be due to the initial conditions of the 
simulations. In this section, we examine the d ifferences be- 
tween o u r simulation s and those performed by INipoti et al.l 
(|2009bl ') . INipoti et"al. ( .2009b. ) construct initial galaxy mod- 
els which have a higher dark matter fraction within the stel- 
lar half mass radius than our models. The correct dark mat- 
ter fraction with the half mas s radii is uncertain for the com - 
pact massive ETGs at high-z . [Sommer- Larsen fc Tof3 (|201G| ) 
show that massive compact high-z galaxies are strongly 
baryon dominated in their inner regions thro ugh cosmolog- 
ical si mulations. The semi-analytic model of IShankar et al] 
(I2OIII ) also predict that the dark matter frac tion of com- 
act m assive ETGs become smaller in high-z. INipoti et alj 
2009bl ) also assume a steeper inner slope, p oc r~^'^ , than 
the Hernquist profile (p oc r~^) for the stellar component. 
Their initial stellar model may preserve the inner structure 
during dry mergers, and may bind the satellite galaxies in 
the inner region. These lead to suppression of significant 
size growth and a lo wer decrease in the velocity dispersion. 
iNipoti et all (|2009bl ) assume more compact satellite galaxies 
than model B of our simulation. This is because they scale 
the satellite galaxies based on the fundamental plane de- 
rived by the SLAGS survey. For model B in our simulations, 
we chose to simply scale down the primary galaxy (model 
A) with the same average density. This is because we as- 
sume that formation epoch of model B is the same as that 
of model A, that is z = 2.0, and that these models are ex- 
pected to have the same average density. This assumption 
may mean that in the dry minor merger simulations, the 
effect of dynamical friction heating may be more effective 
than in our results. 
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6 CONCLUSIONS 

We investigate the effects of sequential dry minor mergers 
on the size evolution of the compact massive ETGs. From 
an analysis of the Millennium Simulation Database, we show 
that such minor mergers are highly likely in the hierarchical 
structure formation. We perform N-body simulations of se- 
quential minor mergers with parabolic and head-on orbits, 
including a dark matter and stellar component, and compare 
the properties of the merger remnants with that from simul- 
taneous minor mergers. We show that the sequential minor 
mergers are the most efficient for size growth and in decreas- 
ing the velocity dispersion of the stellar system. The evolu- 
tion of stellar size and velocity dispersion of the galaxy in the 
sequential minor mergers with compact satellites agrees with 
recent observations. Furthermore, we construct the merger 
histories of candidates of the high-z compact massive ETGs 
using the Millennium Simulation Database, and model the 
size growth of the galaxies via dry minor mergers. We can 
reproduce the mean size growth factor between 2 = 2 and 
2: = if we adopt the most efficient size growth efficiency 
a = 2.7 obtained in the case of sequential minor mergers 
in our simulations. Because these minor mergers are more 
common for the galaxies in more massive FOF haloes (Fig. 
[2|, our results suggest that for ETGs in more massive dark 
matter haloes, the size growth by the dry minor mergers 
is still more effective. However, we note that our numerical 
result is valid for merger histories with typical mass ratios 
between 1/20 and 1/10 and with parabolic and head-on or- 
bits. In addition, our most efficient size growth efficiency is 
likely to an upper limit since it is based on the simulations 
in which the satellites have smaller mass ratios and are on 
radial orbits. In future studies, we will study various cases of 
mass ratios and merger orbits for a more robust prediction. 

In this study, we did not fo cus the "tightness" of the 
local scaling relations of ETGs. iNipoti et al.l (|2009bl ) note 
that it is difficult to explain such an exact relation through 
the mass and size growth of the compact massive ETGs at 
high-z with dry mergers. We will investigate the strength of 
the correlation during the cosmological merger histories in 
our future studies. 
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Figure 1. Loft: Number of our sample galaxies derived by the Millennium Simulation Database (solid lines). Horizontal axis shows the 
bulge mass of the galaxies at 2 = 0. Vertical axis shows the number of galaxies in each mass bin. The mass bins are [1.0 X 10^^ ; 2.0 X 10^^], 
[2.0 X 10" : 5.0 X 10"], and [5.0 X lO" : 2.0 X lO^^] in units of Mq. Fraction of minor merger- dominated galaxies (see text) are also 
plotted (dashed lines with percentages). Right: The same as the left panel, but for the horizontal axis that shows the FOF halo mass of the 
galaxies at 2 = 0. The mass bins are [1.0 X 10^^ . 2.0 x lO^^], [2.0 X lO" : 7.0 X lO"], [7.0 X lO^^ : 2.0 X 10^"'], and [2.0 X lO" : 2.0 X lO^^] 
in units of Mq. 
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Figure 2. Left: Fractional cumulative mass growth by minor mergers with mass ratios M2/M-1 < 1/4 and M2/M1 < 1/10 from 2 = 2.07 
to 2 = for our sample galaxies derived by the Millennium Simulation Database. Here, the fractional cumulative mass growth is the 
fraction of mass increased by the minor mergers from 2 = 2.07 to 2 = to the mass of the bulge at 2 = 0. Horizontal axis shows the bulge 
mass at 2 = 0. Red crosses represent the fractional cumulative mass growth with mass ratios M2/M1 < 1/4. We average these in the same 
mass bins of Fig. [T] and show the average values with black crosses and the solid line. Blue open circles represent that with mass ratios 
M2/M1 < 1/10. We also average these and show them with black open circles and the dashed line. Right : The same as the left panel, but 
for the horizontal axis that shows the FOF halo mass of the galaxies at 2 = 0. 
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Figure 3. Stellar merger histories of our two sample galaxies derived by the Millennium Simulation Database These have never experienced 
major mergers. Left: a central galaxy of a FOF halo which mass is 1.5 x 10^^ Mq. Right: a member galaxy of a FOF halo which mass is 
1.8 X IO^^Mq. Thick solid line is stellar mass of the galaxies as a function of redshift. Bars in eaxii panel show increased mass by minor 
mergers at each event. Dashed and dotted lines show cumulative increased stellar mass from z = 2.07 to z = as minor mergers and 
accretions, respectively. We define a<;cretions as mergers with mass ratios less than 1:20. 
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Figure 4. Density profiles of tlie stellar system and the dark 
matter fialo of a single model A galaxy (Table [TJ simulated as 
a test run to check whether our two-component galaxy model 
is stable. The solid line shows the initial profile {t = 0). The 
dashed line shows the profile after 3.6 Gyr. The concentrated 
ones represent the stellar system. 
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Figure 5. Snapshots of the particles of the three runs onto the x-y plane. Left-hand panels: the simultaneous minor mergers of compa<;t 
satellites, run lAlOBsm. Middle panels: the sequential minor mergers of compact satellites, run lAlOBsq. Right-hand panels: the sequential 
minor mergers of diffuse satellites, run lAlOCsq. 
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Figure 6. The anglc-avcragcd density profiles of the remnant stellar systems of runs. From top left to bottom right: lAlOBsq, lAlOCsq, 
lAlODsq, lAlOBsm, and lAlOCsm. Thiek solid lines are the density profiles of all stars of the merger remnants, dotted lines are of the 
stars that belonged to the primary galaxy (we call primary stars) , and dashed lines are of the stars that belonged to the satellite galaxies 
(we call satellite stars). Thin lines are of the initial galaxy model. The slops of r~^ '^ and r"* *^ are also shown for comparison. Thick arrows 
are half-mass radii of the stellar system of the remnants. Thin arrows are half-mass radii of the stellar system of the initial primary galaxy 
model. 
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Figure 7. Cumulative mass distribution of tlic remnant stellar systems and dark matter particles. From top left to bottom right: lAlOBsq, 
lAlOCsq, lAlODsq, lAlOBsm, and lAlOCsm. Thick solid lines are the cumulative mass distributions of all stars or all dark matter particles 
of the merger remnants, dotted lines are of the primary stars or the primary dark particles. Thin lines are of the initial galaxy model. 
Thick arrows are half-mass radii of the stellar system of the remnants. Thin arrows axe half-mass radii of the stellar system of the initial 
primary galaxy model. 
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Figure 9. Top-left; Linc-of-sight velocity dispersion profiles of 
the remnant stellar systems, lAlOBsm(dottcd), lA10Bsq(thick 
solid), and initial model A (thin solid). Top-rightt: The same 
as the top-left panel, but for lAlOCsm(dotted), lA10Csq(thick 
solid), and initial model A (thin solid). Bottom: anisotropy pa- 
rameter /3 for the merger remnants in the top panels. 



Figure 10. Left-hand panels: initial (top panel) and final (bot- 
tom panel) differential energy distributions of the star particles 
for the simultaneous minor mergers of compact satellites (run 
lAlOBsm). The thick solid, thin solid, and dashed line refer to 
the distributions of total stars, primary stars, and satellite stars, 
respectively. The energy per unit mass E is normalized to the 
internal unit of our simulations. Right-hand panels: same as left- 
hand panels, but for the sequential minor mergers of compact 
satellites (run lAlOBsq). 
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Figure 11. Stellar mass-size relation of the merger remnants and 
the initial galaxy model. Solid triangle is for the initial galaxy 
model A. Empty square is for the run 2A. Solid square is for the 
run lAlOBsm. Empty circle is for the run lAlOBsq. Filled circle 
is for the run lAlOCsm. Empty triangle is for the run lAlOCsq. 
The so lid line shows th e observed stellar mass-size relation for 
ETGs llShen et al.ll2003l) with the dispersion by the dotted lines. 



Figure 13. The average size growth factor (R^ f/R^^i) of our 
sample galaxies derived by the Millennium Simulation Database 
for various size growth efficiencies we assume. The derivation of 
{Re,f / Re,i) a-re described in the text. Horizontal axis shows the 
size growth efficiency a. Empty square is BCGs. Solid triangle is 
normal ETGs. 
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Figure 12. Stellar mass- velocity dispersion relation of the merger 
remnants and the initial galaxy model. Symbols are the same as in 
Figure rm The solid line show s the observed stella r mass- velocity 
dispersion relation for ETGs jNipoti et al.ll2009d) with the one- 
sigma scatter by the dotted lines. 
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Figure 14. The size growth factor R^j/Re-i = {A'lf/Ah)" of 
our sample BCGs (red crosses) and normal ETGs (blue open 
circles) derived by the Millennium Simulation Database, where 
we assume a = 2.7 which is the most efficient case in our results. 
Horizontal axis shows the bulge mass at 2 = 0. Vertical axis shows 

Re,f /Re,i- 
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Figure 15. The same as Fig. 1141 but for tiie horizontal axis that 
shows the FOF halo mass which the galaxies belong to at 2; = 0. 



